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Abstract 



We have searched for two-body charmless decays of B mesons to purely 
hadronic exclusive final states including u> or <j) mesons using data collected 
with the CLEO II detector. With this sample of 6.6 x 10 6 B mesons we 
observe a signal for the u)K + final state, and measure a branching fraction 
of B(B + -> loK + ) = (l.Sto^ ± 0.2) x 10~ 5 . We also observe some evidence 
for the 4>K* final state, and upper limits are given for 22 other decay modes. 
Lhese results provide the opportunity for studies of theoretical models and 
physical parameters. 
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In the last several years, the study of charmless non-leptonic decays of B mesons has 
attracted a lot of attention, primarily because of the importance of these processes in under- 
standing the phenomenon of CP violation. This interest is expected to continue as several 
new experimental facilities specifically built for B meson studies begin operating within a 
few years. Purely hadronic decays of B mesons are understood to proceed mainly through 
the weak decay of a h quark to a lighter quark, while the light quark bound in the B meson 
remains a spectator, as shown by the Feynman diagrams in Figure |]. The decay amplitude 
for "tree- level" b — > u transitions (Figure |I](a) and (b)) is much smaller than the one for dom- 
inant b — > c transitions due to the ratio of Cabibbo-Kobayashi-Maskawa [|l|] matrix elements 
Vub/Vcb ~ 0.1. Transitions to s and d quarks are effective flavor- changing neutral currents 
proceeding mainly by one-loop "penguin" amplitudes, and are also suppressed. Examples 
are shown in Figure |I] (c) and (d). 
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FIG. 1. (a) and (b) tree-level spectator and (c) and (d) penguin diagrams for some of the decay 
modes investigated. 

The strong interaction between particles in the final state makes theoretical predictions 
difficult. The use of effective hamiltonians, often with factorization assumptions PHiOl, 
has led to a number of these predictions, and the experimental sensitivity has now become 



sufficient to allow us to begin to test the correctness of the underlying assumptions ||TTHT3 

In this letter, we describe searches for I?-meson decays to exclusive final states that 
include an u or <ft meson and one other low-mass charmless meson. Some decays to final 
states with a are of particular interest because they are dominated by penguin amplitudes, 
and receive no contribution from tree- level amplitudes (see Figure [I]), while others, such 
as B + — > 07r + , receive no contribution from penguin or tree amplitudes and only proceed 
through higher-order diagrams. 
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The results presented here are based on data collected with the CLEO II detector |TJ| at 
the Cornell Electron Storage Ring (CESR). The data sample corresponds to an integrated 
luminosity of 3.11 fb" 1 for the reaction e + e~ — > T(AS) — > BB, which in turn corresponds 
to 3.3 x 10 6 BB pairs. To study background from continuum processes, we also collected 
1.61 fb™ 1 of data at a center-of-mass energy below the threshold for BB production. 

The final states of the decays under study are reconstructed by combining detected 
photons and charged pions and kaons. The u and <fi mesons are identified via the decay 
modes uj — > 7r + 7r _ 7r° and <fi — > K + K~ , respectively. The detector elements most important 
for the analyses presented here are the tracking system, which consists of 67 concentric drift 
chamber layers, and the high-resolution electromagnetic calorimeter, made of 7800 CsI(Tl) 
crystals. 

Reconstructed charged tracks are required to pass quality cuts based on their track fit 
residuals and impact parameter. The specific ionization (dE/dx) measured in the drift layers 
is used to distinguish kaons from pions. Expressed as the number of standard deviations 
from the expected value, Si(i = tt,K), it is required to satisfy < 3.0. Photons are defined 
as isolated showers, not matched to any charged tracks, with a lateral shape consistent with 
that of photons, and with a measured energy of at least 30 (50) MeV in the calorimeter 
region | cos^l < 0.71(> 0.71), where 6 is the polar angle. 

Pairs of photons (charged pions) are used to reconstruct 7r°'s and 77's (K 0> s). The mo- 
mentum of the pair is obtained with a kinematic fit of the decay particle momenta with 
the meson mass constrained to its nominal value. To reduce combinatoric background, we 
reject very asymmetric ir° and 77 decays by requiring that the rest frame angle 9* between the 
direction of the meson and the direction of the photons satisfy | cos#*| < 0.97, and require 
that the momentum of charged tracks and photon pairs be greater than 100 MeV/c. 

The primary means of identification of B meson candidates is through their measured 
mass and energy. The quantity AE is defined as AE = E\ + E% — Eb, where E\ and E2 are 
the energies of the two daughter particles of the B and Ef, is the beam energy. The beam- 



constrained mass of the candidate is defined as M = J El — |p| 2 , where p is the measured 
momentum of the candidate. We use the beam energy instead of the measured energy of the 
B candidate to improve the mass resolution by about one order of magnitude. 

The large background from continuum quark-antiquark (qq) production can be reduced 
with event shape cuts. Because B mesons are produced almost at rest, the decay products 
of the BB pair tend to be isotropically distributed, while particles from qq production have 
a more jet-like distribution. The angle 9t between the thrust axis of the charged particles 
and photons forming the candidate B and the thrust axis of the remainder of the event is 
required to satisfy | cos#t| < 0.9. Continuum background is strongly peaked near 1.0 and 
signal is approximately flat for this quantity. We also form a Fisher discriminant (J 7 ) fllT| 
with the momentum scalar sum of charged particles and photons in nine cones of increasing 
polar angle around the thrust axis of the candidate and the angles of the thrust axis of the 
candidate and p with respect to the beam axis. 

The specific final states investigated are identified via the reconstructed invariant masses 
of the B daughter resonances. For final states with a pseudoscalar meson, and for the 
secondary decay 77' — > pj, further separation of signal events from combinatoric background 
is obtained through the use of the defined angular helicity state of the <fi, uj, or p. The 
observable TC is the cosine of the angle between the flight direction of the vector meson 
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and the daughter decay direction (normal to the decay plane for the u), boosted to the 
meson's rest frame. For the final states uK* + and up + , the 7r° from K* + or p + decay defines 
the daughter direction. In this case we require Ti < 0.5 to reduce the large combinatoric 
background from soft 7r°'s. Since the distribution of 7i is not known for these vector- vector 
final states we assume the worst case {Ti 2 ) when computing the efficiency 

Signal event yields for each mode are obtained with unbinned multi-variable maximum 
likelihood fits. We also performed event counting analyses that applied tight constraints on 
all variables described above. Results for the latter are consistent with the ones presented 
below. 

For N input events and p input variables, the likelihood is defined as 

N P 

C = e -(N 3+ N B) | I{ Vs 1J Vs { fir .... fm , Xij) + 
i=l j=l 

V 

NBllVB^giji-^nj^ij)}, 

i=i 

where Vs %j and Vb v are the probabilities for event i to be signal and continuum background 
for variable Xij, respectively. The probabilities are also a function of the parameters / and 
g used to describe the signal and background shapes for each variable. The number of 
parameters required varies depending on the input variable. The variables used are AE, M, 
T, resonance masses, and TC as appropriate. For pairs of final states differentiated only by 
the identity of a single charged pion or kaon, we also use Si for that track and fit both modes 
simultaneously. N$ and Nb, the free parameters of the fit, are the number of signal and 
continuum background events in the fitted sample, respectively. We verified that background 
from other B decay modes is small for all channels investigated and did not require inclusion 
in the fit. Correlations between input variables were found to be negligible, except between 
the invariant masses of a parent resonance and its daughter, which the likelihood function 
takes into account. 

For each decay mode investigated, the signal probability distribution functions (PDFs) 
for the input variables are determined with fits to Monte Carlo event samples generated 
with a GEANT (IT5J based simulation of the CLEO detector response. The parameters of 



the background PDFs are determined with similar fits to a sideband region of data defined 
by \AE\ < 0.2 GeV and 5.2 < M < 5.27 GeV/c 2 . The data samples collected on and below 
the T(4S) resonance are used. The signal shapes used are Gaussian, double Gaussian, and 
Breit-Wigner as appropriate for AE and mass peaks. For background, resonance masses are 
fit to the sum of a smooth polynomial and the signal shape, to account for the component 
of real resonance as well as the combinatoric background. For AE and M background we 
use a first-degree polynomial and the empirical shape f(z) oc My/1 — z 2 exp (— £(1 — z 2 )), 
where z = M/E b and £ is a parameter to be fit, respectively. Finally, for T , Sk, and S n , we 
use bifurcated Gaussians for both signal and background. 

Sideband regions for each input variable are included in the likelihood fit. The number of 
events input to the fit varies from 70 to ~ 12000, depending on the final state. Table | [IE| 



gives the results for each mode investigated. The final state uh + represents the sum of the 
ujK + and uotx + states (h + = K + or tt + ). Shown are the signal event yield, the efficiency, 
the product of the efficiency and the relevant branching fractions of particles in the final 



6 



state, and the branching fraction for each mode, given as a central value with statistical and 
systematic error, or as a 90% confidence level upper limit. The one standard deviation (er) 
statistical error is determined by finding the values where the quantity x 2 — — 2 ln(£/£ max ), 
where £ max is the point of maximum likelihood, changes by one unit. 

Systematic errors are separated into two major components. The first is systematic errors 
in the PDFs, which are determined with a Monte Carlo variation of the PDF parameters 
within their Gaussian uncertainty, taking into account correlations between parameters. 
The final likelihood function is the average of the likelihood functions for all variations. The 
second component is systematic errors associated with event selection and efficiency factors. 
For cases where we determine a branching fraction central value, the final systematic error 
is the quadrature sum of the two components. For upper limits, the likelihood function 
including systematic variations of the PDFs is integrated to find the value that corresponds 
to 90% of the total area. The efficiency is reduced by one standard deviation of its systematic 
error when calculating the final upper limit. 

For final states which we detect in multiple secondary channels, we sum the value of \ 2 
as a function of the branching fraction and extract the final branching fraction or upper 
limit from the combined distribution. Table [TT] shows the final results, as well as previously 
published theoretical estimates. 

We find a significant signal for B + — ► ujK + and measure the branching fraction B(B + — > 
uK+) = (l.Btol ± 0-2) x 10~ 5 , where the first error is statistical and the second system- 
atic. We also find a signal for B + — > ujh + , with a branching fraction of B{B + — * uh + ) = 
(2.5lo!7 =t 0.3) x 10~ 5 . The significance for these signals is 3.9a for B + — > ujK + and 5.5a 
for B + — > ojh + . We also find some evidence for the sum of the modes B + — > <j)K* + and 
B° — > (f)K*°, with a significance of 2.9a. It is sensible to combine these modes since their 
decay rate is expected to be dominated by identical penguin amplitude contributions, except 
for different spectator quarks. The quoted significances include both statistical and system- 
atic errors. If we interpret the observed <pK* event yield as a signal, we obtain a branching 
fraction of B(B -> (j)K*) = (l.l±g:| ± 0.2) x 10~ 5 . Figure | shows the likelihood functions 
for these modes. Figure [3] shows the projection along the M axis, with clear peaks at the B 
meson mass. 

We also set lower limits on the branching fractions for B + — > ujK + and B + —>■ uh + , which 
could have interesting theoretical implications |T7-2"I|. We find B(B + — ► ujK + ) > 8.4 x 10~ 6 
and B(B + — > ujh + ) > 1.6 x 10~ 5 at the 90% confidence level. The latter limit would imply 
that the parameter £ used in references |T3 and ETJ is restricted to the region £ > 0.62 



and £ > 0.53, respectively. However, based on reference [H| our measurement of B(B + — > 
cf)K + ) < 0.5 x 10 -5 implies that £ < 0.27 at the 90% confidence level. Although there is 
still considerable uncertainty in the theoretical model parameters, these limits illustrate the 
difficulty in accounting for all our current results with a single phenomenological parameter. 

We thank A. Ali, H. Lipkin, J. Rosner, H.-Y. Cheng, and S. Oh for useful discussions. We 
gratefully acknowledge the effort of the CESR staff in providing us with excellent luminosity 
and running conditions. This work was supported by the National Science Foundation, the 
U.S. Department of Energy, Research Corporation, the Natural Sciences and Engineering 
Research Council of Canada, the A. P. Sloan Foundation, and the Swiss National Science 
Foundation. 
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TABLE I. Measurement results. Columns list the final states (with secondary decay modes as 
subscripts), event yield from the fit, reconstruction efficiency e, total efficiency including secondary 
branching fractions B s , and the resulting B decay branching fraction B. 



Final state 


Yield(events) 


e(%) 


eB s {%) 


B(10~ 5 ) 


luK + 


12 2 +5 A 

- LZ " Z '— 4.5 


28 


25.1 


1 5+S-I ± 

"—0.6 u 


ujK° 


1.5 


15 


4.4 


< 5.7 




q o+ 5 - 3 
y - z -4.3 


9Q 






uih + 

UJ 1 1 


21 4+ 6 - 5 


99 


95 5 




LJ1T 


2 4+?' 9 


24 


20.9 


< 1.4 




o i +1 - 9 

u - 1 -o.i 


16 

-L U 


2.4 


<r 6 4 


Ml/ 


5 i + H 


16 


4.2 


< 9.2 




—0.0 


24 


8.5 


< 2.0 




+ 2'£ 


15 


3.2 


< 2.8 


UjK*t n 


J --- L — 1.1 


7 


2.0 


< 12.9 




4 5+H 

^•^—2.8 


16 


3.2 


< 10.9 




2 

i,J --2.1 


22 


13.1 


< 2.3 


LO p^ 


2 5+M 

^•"—2.5 


8 


6.8 


< 6.1 




o o + H 


24 


21.1 


< 1.1 


COLO 


u - J -0.3 


15 


11.9 


< 1.9 




n n+ - 8 


47 


9^ 1 




4>K° 


1 q+2.0 


32 


5 3 


< 3 1 




n o +0 - 9 
u - u -o.o 




94 n 






o o+S-S 

u u — 0.0 


31 


15.1 


< 0.5 




o.ol8:S 


26 


2.2 


< 3.5 


Mn 


2.71U 


30 


4.4 


< 6.3 




n o +0 - 6 

u - u -o.o 


39 


7.5 


< 1.3 


<Pmn 


n o +0 - 5 

u - u -o.o 


24 


2.7 


< 2.9 




2-6l|l 


26 


4.4 


< 5.6 




i.7±?:S 


29 


3.4 


< 5.3 




3-2l|? 


39 


12.7 


< 2.2 




o o +L9 


18 


1.0 


< 8.0 


<f>p+ 


0-Olao 


34 


16.7 


< 1.6 


<pp° 


o-si^i 


41 


20.0 


< 1.3 






23 


10.2 


< 2.1 




0.4l 1 


40 


9.7 


< 1.2 
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TABLE II. Combined results and expectations from theoretical models. 



Decay mode 


0(io~ 5 ) 


Theory B (10~ 5 ) 


References 


B+ 


-> ojK + 


1.5^ ±0.2 


0.1 - 0.7 


§§ 


SIC] 


B°- 


^uK° 


< 5.7 


0.1 - 0.4 


11, 


10 


i 




B+ 


— ► uir + 


< 2.3 


0.1 - 0.7 


II 


ETC] 


B+ 


-► uh+ 


2.5lg:?±0.3 


- 










B°- 




< 1.4 


0.01 - 1.2 


II 


I 


i 




B°- 


-* urj' 


< 6.0 


0.3 - 1.7 


IP 

1 II 


I 






B°- 


-* UJT] 


< 1.2 


0.1 - 0.5 


§i 


1 






B+ 


-» ojK* + 


< 8.7 


0.04 - 1.5 


B§ 


1 






B°- 




< 2.3 


0.2 - 0.8 


| llll 








B+ 


— > 


< 6.1 


1.0-2.5 


§§ 


1 






B°- 


- V 


< 1.1 


0.04 


1 








B°- 




< 1.9 


0.04 - 0.3 


B§ 








B+ 




< 0.5 


0.07- 1.6 


Hi 


1 




S0 


B°- 




< 3.1 


0.07- 1.3 


§1 


M 


,1C] 


B+ 


— * (f)ir + 


< 0.5 


« 0.1 


m 


|, 


m 


B°- 




< 0.5 


« 0.1 


M 


I 


i 




B°- 


-> 07/' 


< 3.1 


« 0.1 


EE 


I 






B°- 


-» 07/ 


< 0.9 


« 0.1 


m 


1 


i 




B+ 


-> 0K*+ 


< 4.1 


0.02 - 3.1 


m 


10 


1 


B°- 




< 2.1 


0.02 - 3.1 


ii 


10 




B - 


► c^RT* 


< 2.2 


0.02 - 3.1 


@s 


m 




B+ 




< 1.6 


« 0.1 


ii 


i 






B°- 




< 1.3 


« 0.1 


ii 








B°- 


-> 0a; 


< 2.1 


« 0.1 


ii 








B°- 




< 1.2 


none 
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FIG. 2. (a) Likelihood function contours for B + — ► coh + ; (b) The function 
-21n£/£ max = x 2 - X^n for 5 - <\>K\ 
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5.30 



B mass (GeV / c ) 



B mass (GeV / c ) 



FIG. 3. Projection onto the variable M for (a) B + — ► loK + (shaded) and B + — > w7r + (open) and 
(b) -B — > </>iiC*. The solid line shows the result of the likelihood fit, scaled to take into account the 
cuts applied to variables not shown. The dashed line shows the background component, and in (a) 
the dotted line shows the B + — > ujK + fit only. 
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